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Abstract 

Predictions for the radiative return with the muon pair and pion pair final state 
from KKMC and PHOKHARA Monte Carlo programs are compared and discussed. 
The case of muon pairs is well understood, especially of the initial state radiation 
(ISR), where three different second order calculations agree very well. The case of 
the final state radiation (FSR) requires more tests. Matrix element in KKMC of the 
EEX type with the incomplete second order NLL corrections is not good enough 
for the radiative return at Q 2 < lGeV with the precision requirement better than 
1%. A method of extending the superior CEEX-type matrix element in KKMC to 
the pion pair final state is described. 
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1 Introduction 



The aim of this contribution is to compare whatever the best we have at hand for evalua- 
tion of the initial state radiation (ISR) effect in the process e _ e + — > fi~ using KKMC 
p] |2j and PHOKHARA EJ IE] Monte Carlo programs. The above investigation will be 
partly extended to the process e~e + — > 7r~7r + 7. 
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Figure 1: Mu-pair mass (square) spectrum in case of ISR only. ^Js = 1.01942GeV. 



2 ISR in muon pair production 

Both KKMC and PHOKHARA programs are full scale MC event generators, which can 
provide for any experimentally observable distribution. We concentrate however, on the 
distribution of the squared mass spectrum Q 2 = s' of the muon pair, because this distri- 
bution is relevant for the radiative return measurements of R(s), and also because this 
particular distribution we may compare with the classical semi-analytical calculations. 
Here we shall also exploit the analytical formulas of ref. |Hj (see also [Zj), which imple- 
ment analytical second order ISR calculation of ref. |Hj and third order leading-logarithmic 
(LL) ISR calculation of refs. The ISR formula of ref. jE] is provided by the KKsem 

facility of KKMC. In the actual KKsem implementation we use version of the formula 
where numerically negligible (at least at LEP energies, see ref. jSj) second order NNLL 
terms are neglected. 



1 



It is important to stress from the very beginning that authors of PHOKHARA and 
KKMC use different terminology to describe Born level matrix element and higher order 
matrix element. I shall not try to unify terminology or fully explain the differences, 
referring the reader to original works, like refs. P3I2] and jSlSIS]- Let me explain only very 
briefly the main differences. The KKMC authors define Born as e + e~ — > ff without any 
photon emission and the radiative return is necessarily the first order process with respect 
to such a Born level. The leading-logarithmic (LL) corrections are of order a n L n , where 
n = 1,2, 3.. .00 is the standard perturbative order, while mass logarithm L = \\i{s/m 2 e ) 
is coming either from the virtual photon correction or the phase space integration over 
the real photon angle down to zero value. The NLL and NNLL corrections are of order 
a n L n ~ x and a n L n ~ 2 correspondingly. Concerning mass terms, they are routinely neglected 
in KKMC for the electron (except those which integrate up to a finite correction) while 
an effort is made to keep all of them for the final state fermions, at least at the Born and 
the first order level. KKMC implements several variants of the QED matrix elements, 
which feature different level of higher order and mass term truncation. PHOKHARA 
authors employ the leading-order (LO) as a name for the process in which one (and only 
one) photon is emitted in the final state. They name as the next-to-leading-order (NLO) 
their matrix element with the one-loop corrections and the second real photon. This 
terminology may seem more adequate to discuss radiative return. However, when trying 
to match the two terminologies one has to pay attentions to the available phase space 
of the first real photon. Depending on whether the minimum emission angle is imposed 
or not, one gets full factor L or not, even at the LO. This affects strongly the relative 
magnitude of higher order corrections with respect to LO or LL. In this study we generally 
exclude from the considerations "non-photonic" corrections due to emission of additional 
lepton pair and vacuum polarization. 

In Fig^ we compare results of KKMC and of PHOKHARA using the best available 
ISR matrix element in both programs at ^/s = 1.01942GeV. In KKMC we use second order 
matrix element with coherent exclusive exponentiation (CEEX) described in refs. |2*| lllj. 
The second order CEEX matrix element has complete next-to-leading-logarithmic (NLL) 
contributions 1 and complete next-next-to-leading-logarithmic (NNLL) contributions. The 
magnitude of NLL and NNLL corrections was also examined in a separate studies, see 
contribution of S. Yost in these proceedings. KKMC includes most of the third order LL 
contributions by the virtue of exponentiation 2 . On the other hand, PHOKHARA imple- 
ments complete second order ISR, including complete NLL and NNLL corrections (i.e. 
singular corrections proportional to ^m 2 e and ^m^, which integrate to finite corrections of 
order -, in the limit m e — > 0). PHOKHARA does not resum (exponentiate) soft photon 
contributions to infinite order. It is worth to stress that the two MC calculation, KKMC 
and PHOKHARA, and semianalytical formula of KKsem of refs. [UJ |H] represent set of 
three completely independent second order (using terminology of KKMC) calculation of 
the ISR in every aspect of calculating QED matrix element and integrating the phase 

1 For unpolarized beams, see discussion in ref. related to eq. (128) therein. 

2 It is also known that exponentiation of the YFS type sums up quite substantial part of third order 
LL, see refs. [TO] 
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space. 

The main comparison of the ISR calculations is shown in FigHt, where the distribution 
da/dQ 2 from KKMC and KKsem agree very well, within 0.2%, except very low Q 2 where 
they diverge by about 0.3% 3 , while Fig^, shows certain addition cross-check. The reason 
for this discrepancy is not clear. Neglected NNLL in KKsem are a viable candidate, but 
to confirm this hypothesis one would need more tests. In the same plot we see that 
PHOKHARA agrees well with KKsem at low Q 2 (aligning with KKsem) and differs by 
about 0.25% in the central region (we need higher statistics from PHOKHARA to confirm 
this number) from both KKMC and KKSEM and drops sharply at soft limit, high Q 2 , 
because of lack of soft photon resummation. In order to understand quantitatively the 
effect of lack of exponentiation in PHOKHARA we compare its result in FigEh with a 
variant of KKsem in which we switch off exponentiation, i.e. all terms beyond second 
order are truncated. The smooth curve in FigCb representing result of this truncation 
fits very well PHOKHARA result. In particular, looking into this result, one may think 
that the deviation of PHOKHARA by 0.25% in the central region is related to its neglect 
of the third order LL. This conjecture needs more test to be confirmed. 
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Figure 2: Muon pair mass spectrum from KKMC and KKsem. 



We summarize on the results of Fig^that KKMC with the second order CEEX matrix 
element, PHOKHARA with its second order matrix element and KKsem implementing 
second order analytical calculation agree very well, within the expected range and the 
pattern of the discrepancies seems to be understood. 

In KKMC there is another more primitive QED matrix element denoted as EEX, 
see ref. [2] for its full description, which follows closely the classical Yennie-Frautschi- 

3 Note that similar comparison of KKMC and KKsem was done in ref. [2] for LEP energies. At the 
present lower energy y/s — 1.01942GeV subleading terms are, however, more important. 
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Suura (YFS) exponentiation scheme and its implementation is limited to first order plus 
second order LL. In the second order EEX matrix element (contrary to CEEX) the NLL 
corrections are incomplete. (On the other hand EEX third order LL is complete, while 
in CEEX it is incomplete.). For technical and historical reasons, see discussion below, 
EEX type matrix element is used for the production of low energy hadronic final states, 
for example for pion pair. It is, therefore, important to check how good it is compared to 
KKMC with more complete coherent exclusive exponentiation (CEEX) matrix element. 
This is done for the muon final state in Figs. El (CEEX is not yet available for re- 
pairs). In Fig. we see results from KKMC CEEX and several variants of EEX. We 
are actually plotting da/dQ 2 , dividing all results by KKsem of ref. jSj, the same as in 
previous FigO The curves marked EEX72 represent exponentiated EEX matrix element 
based on complete first order, while EEX73 and EEX74 include also complete second 
and third LL, while CEEX203 4 is the same as in Fig. ^ As we see, at low Q 2 , that is 
for the hard photon emission, results of EEX matrix elements depart from other more 
complete results by up to 3%! In the p region it is different from the KKsem, CEEX 
KKMC by about 1%. The above result is also consistent with what we have seen in 
Fig. d EEX is therefore not well suited for the use in the high precision measurements 
of R(s) using radiative return below Q 2 =lGeV. This result is not very much surprising, 
as EEX of KKMC has incomplete second order NLL. The observed effect at the low Q 2 
is a little bit bigger then what we expected. We have therefore done certain additional 
tests. We have split EEX results into three components, fa, i = 0, 1, 2, compared each 
of them with analytical result of table I of ref. j2], also at yfs = lOGeV. We do not show 
results of these tests here, but the overall pattern of discrepancies seems to be consistent 
with NLL class of corrections. This additional test indicates also that the main source of 
the problem is an approximate double real emission matrix element in EEX and not the 
incomplete virtual corrections. In particular we have included in these tests the complete 
NLL contribution in fa and fa. This did not help! The whole discrepancy seems to result 
from the use of the LL-approximate matrix element for the double real photon emission in 
EEX. The above observation is consistent with the older tests in ref. j2] at LEP energies. 

2.1 Muon pair, ISR+FSR 

Let us not include FSR in the game, again for the muon pair final state. In fact, at low 
Q 2 the rate of muon pair in radiative return is higher than of tt pairs, hence da/dQ 2 of 
muons can be used as a reference distribution for measuring R(s). It is therefore worth to 
test FSR in KKMC and to check result for da/dQ 2 once again. In Fig. we show result 
from KKMC for second order CEEX matrix element in which we include ISR, FSR and 
its interference. We compare MC result with the semianalytical result of KKsem in which 
with the same ISR radiative function of ref. jB] . The FSR distribution of ref. [2] features 
incomplete NLL in KKsem, so it is definitely inferior with respect to ISR counterpart - 
the complete list of the FSR radiative corrections in KKsem can be found in Table II in 
ref. j2] • This above semianalytical formula also misses the interference of ISR and FSR, 

4 Indices 203, 74 etc. follow numbering of MC weights in ref. [I]. 
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Figure 3: Comparison of PHOKHARA NLO and KKMC with EEX matrix element. ISR only. 



which in first order is zero in the inclusive da/dQ 2 so this omission does not harm. In 
Fig. 12^. we see the ratio of the corresponding results from KKMC and KKsem. (NB. 
PHOKHARA is able to provide result with FSR for muon pairs in the LO, and it would 
be interesting to include it in the comparison.) This result is rather preliminary and has 
to be checked. In any case, the agreement better than 1% found all over the dQ 2 range is 
quite satisfactory as a starting point for further investigation 5 

5 In one bin we see trace of large weight fluctuation which is probably due to rounding errors. This 
result was obtained using weighted events. For the MC run with weight-one events this effect would 
disappear. Such numerical instabilities need further investigation. 
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3 ISR for tt + 7t pair production 

Let us now switch to low Q 2 7r-pair state produced at the radiative return process. In Fig. El 
we compare KKMC with the EEX style matrix element on one hand with PHOKHARA 
second order (marked as PHOKHARA2) on the other hand. The EEX matrix element 
is the default one in KKMC, with first order exponentiation, the completed second and 
third order LL (EEX74). We limit ourselves to ISR only. The results of Figs. E^-b are 
obtained without any cutoffs. In Fig. Ek we show the actual distributions, including also 
the distribution for the muon pairs. We see that for Q 2 < 0.33 the muon-pair cross section 
is bigger than that of 7r-pair. The ratio PHOKHARA/KKMC is not so well understood as 
the the analogous results for the muon pair shown in the previous section. The discrepancy 
at high Q 2 we attribute to lack of exponentiation in PHOKHARA 6 while another larger 
discrepancy at low Q 2 is most likely due to incompleteness of second order NLL in EEX 
matrix element on KKMC, and it corresponds to deviation which was already seen in 
Fig. |2K- In Fig. Efc-d we show the analogous results for relatively mild cut on photon 
momentum, where photon is defined as a "missing four momentum" calculated knowing 
pion momenta and beam momenta. We ask for the momentum of such a "collective 
unseen photon" to be directed below 15° from the beam and to have at least lOMeV 
of energy. For each ir we require that it is situated in wide angles, eg. separated by 
more than 40° from each beam. Results in Fig. Efc-d look quite similar, except that the 
discrepancy between PHOKHARA and EEX KKMC is bigger (we need better statistics 
from PHOKHARA to see it more clearly). This can be attributed to the fact that the 
leading logarithm L due to real emission is diminished by the cut on the photon angle 
with respect to beams. 

3.1 How to extend CEEX ISR to hadronic final states? 

In the following we show that the superior CEEX ISR matrix element can be extended to 
hadronic final states at low Q 2 , like pion pair. This can be done provided we have some 
decent modeling of the hadronic final state in terms of the corresponding formfactor. In 
view of its practical importance, let us elaborate on this point. 

6 This does not hinder practical applications of PHOKHARA for radiative return measurements, which 
concentrate at lower Q 2 . 
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In CEEX Born amplitude for ee — > /i/z is defined as a four-spinor tensor 



ie 2 



B=i,Z B=j,Z 

(Gf^iba] = v{p b , Xb)Gf tfl u(p a , A ), (C^)^] = w(p c , A c )G^<u(p d) A d ) 



(1) 



A=± A=± 

fJbV 



11% {X) 



X 2 - M B 2 + iT B X 2 /M B 1 

and it enters as a basic building block in every spin amplitude in the CEEX scheme, with 
arbitrary number of photons. See eq. (43) in ref. [2] for notation. The above Born is 
calculated using Chisholm identity and replaced with the bi-spinor objects of the Kleiss- 
Stirling method. 

In case of hadronic final state the structure of the Born amplitude is 

% a] (X)J,(X,q t ), J,(X, gi )X» = 0, (2) 
where qi are momenta of the final state hadrons, X = ^2 an d 

% a] (X)=ie 2 Yl H * (OH n *P0 = E ®m(*)- (3) 

B=y,Z B=y,Z 

— * 

In the rest frame of X one has = (0, J) and we may split J into difference of the 
two massless four- vectors = J+ — Jt. In the arbitrary reference frame the above 
prescription extends as follows 7 

1 



J£ = ^=(V^PX»±VX>J»). (4) 

Each of the two corresponding components in 03 can be expressed in terms of the of the 
Kleiss-Stirling bi-spinors s±(pi,p 2 ), see eqs. (A4-A6) in ref. This can be done using 
completeness relation for v(pb, A&) y u J± u (.Pa, A a ) taking advantage of J\ = 0. 

Note that the above CEEX implementation requires that we parametrize the produc- 
tion amplitude of the each final hadronic state one by one in terms of the formfactors in 
a completely exclusive manner. However, this is necessary anyway for good phenomeno- 
logical description of these often resonant low energy hadronic states. We conclude that 
CEEX can be used for ISR for low energy hadron production. The question is only how 
much programming it will be and who will do it. 

7 The author would like to thank J. Kuchn for suggesting him this solution. 
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3.2 Conclusions 



The case of radiative return with the muon pair final state is well understood, especially 
for the ISR where three different second order calculations agree very well. The case with 
FSR requires more tests. Matrix element EEX of KKMC with the incomplete second 
order NLL is not good in the radiative return at Q 2 < lGeV for precision requirement 
better than 1%. Method of porting CEEX matrix element of KKMC to pion pair final 
state is outlined. 
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